Spin splitting and precession in quantum dots with spin-orbit coupling: the role of 

spatial deformation 
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Extending a previous work on spin precession in GaAs/AlGaAs quantum dots with spin-orbit 
coupling, we study the role of deformation in the external confinement. Small elliptical deformations 
are enough to alter the precessional characteristics at low magnetic fields. We obtain approximate 
expressions for the modified g factor including weak Rashba and Dresselhaus spin-orbit terms. For 
more intense couplings numerical calculations are performed. We also study the influence of the 
magnetic field orientation on the spin splitting and the related anisotropy of the g factor. Using 
realistic spin-orbit strengths our model calculations can reproduce the experimental spin-splittings 
reported by Hanson et al. (cond-mat/0303f39 ) for a one-electron dot. For dots containing more 
electrons, Coulomb interaction effects are estimated within the local-spin-density approximation, 
showing that many features of the non-iteracting system are qualitatively preserved. 

PACS numbers: 73.21. La, 73.21.-b 



I. INTRODUCTION 

In the last years the study of spin-related phenomena 
has become one of the most active research branches in 
semiconductor physics. The present advances in spin- 
based electronics^ and the hope for better devices, with 
enhanced performance with respect to the conventional 
charge-based ones, encourage this research. Two physi- 
cal mechanisms underlie the operation of most spintronic 
devices: a) the spin-spin interaction, present in ferromag- 
netic materials and in diluted magnetic semiconductors; 
and b) the electron spin-orbit (SO) coupling stemming 
from relativistic corrections to the semiconductor Hamil- 
tonian. It should also be mentioned that, as shown re- 
cently by Ciorga et al., another possibility of spin con- 
trol involves the use of external magnetic fields to induce 
changes in the spin structure of a quantum dot«£ These 
spin modifications affect the passage of currents through 
the system, originating the spin blockade effect. A con- 
spicuous example of device exploiting the SO coupling is 
the spin transistor, first proposed by Datta and Das^ In 
this system the spin rotation induced by an adjustable 
Rashba coupling is used to manipulate the current. 

In a recent work^ we studied the spin precession of 
quantum dots with SO coupling under the action of a 
vertical magnetic field of modulus B. It was shown that 
the SO coupling modifies the precessional frequency from 
the Larmor expression Hlvl = \9*\hbB, where g* is the 
bulk effective g factor and /ib is the Bohr magneton, to 
a different value depending on the dot quantum state. 
Namely, the modified precessional energy equals the gap 
between spin up and down states for the active level, the 
so called spin-flip gap A s /. 

Purely circular dots are characterized by discontinuous 
jumps in angular momentum with the number of elec- 
trons N and the magnetic field, with a similar behavior 
for the precessional frequency. An interesting prediction 
of Ref. was that for some values of N a finite A s / per- 
sists even at B = 0, i.e., a constant offset to the above 



Larmor formula. It is our aim in this work to extend 
those investigations by including deformation in the ex- 
ternal confinement, as well as a more general treatment 
of the SO coupling, considering Rashba and Dresselhaus 
contributions on an equal footing. We shall show that 
small elliptical deformations are enough to sizeably al- 
ter the precessional frequency, yielding a deformation- 
dependent g factor and washing out the low B offsets of 
purely circular dots. Anisotropy effects in the g factor 
will also be studied by allowing for a tilted orientation of 
the magnetic field vector with respect to the dot plane. 

Spin dynamics in semiconductor nanostructures can be 
experimentally monitored with optical techniques. In- 
deed, a time-delayed laser interacting with a precessing 
spin experiences the Faraday rotation of its polarization. 
Measuring the rotation angle for different delays allows 
to map the spin orientation and thus observe in detail 
the dynamics. This technique has been applied to bulk 
semiconductors (see Ref. for a recent review) and, also, 
to CdSe excitonic quantum dots in Ref. 0. Alternative 
methods to gather information on the g factor in quan- 
tum dots normally use measurements of the resonant tun- 
neling currents through the system that permit the de- 
termination of the spin splittings and, therefore, deduce 
the effective g value i£ 

Electron spin in quantum dots is much more stable 
than in bulk semiconductors, due to the suppression of 
spin flip decoherence mechanisms^ Spin relaxation is pre- 
dicted to occur on a time scale of 1 ms for B = 1 T. Ac- 
cordingly, in this work we shall neglect spin relaxation, 
focussing on the much faster spin precession in quantum 
dots. The spin splittings will be compared with those 
measured in Ref. |g, showing that realistic values of the 
SO strengths can indeed reproduce the observed behav- 
ior. The paper is organized as follows. Section II presents 
the analytical model for low SO intensities. In Sec. Ill we 
discuss the numerical results for a variety of situations; 
namely, arbitrary SO strengths (A), tilted magnetic fields 
(B), one-electron dots (C) and treating Coulomb inter- 
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action effects (D) within the local-spin-density approxi- 
mation (LSDA). Finally, the conclusions are presented in 
Sec. IV. 



II. THE MODEL 

A. The noninteracting Hamiltonian 

Our model of a single quantum dot consists in TV elec- 
trons of effective mass m* whose motion is restricted to 
the xy plane where an electrostatic potential V ex t(r) in- 
duces the confinement. We assume a GaAs host semicon- 
ductor, for which m* = 0.067m e . To allow for elliptically 
deformed shapes we consider an anisotropic parabola, 
i.e., 



V«t(r) = -m*(iv 2 x x 2 +u 2 y y 2 ) 



(1) 



Neglecting for the moment Coulomb interactions between 
electrons we treat the Hamiltonian for independent par- 
ticles Hip = X^iliM*)- The single-electron Hamilto- 
nian (h) contains the kinetic/confinement energy (h ), 
the Rashba (hu) and Dresselhaus (ho) SO terms and 
the Zeeman energy (hz)] 



h = h + h R + h D + h z ■ 
The explicit expressions of ho and hz read 



h = 



2m* 
1 



+ Vext(x,y) 



hz = ~g* Hb(B x (J x + B y Oy + B z a z ) 



(2) 

(3) 
(4) 



where P = —ihV + - A represents the canonical momen- 
tum depending on the vector potential A = B z /2(—y, x) 
and the cr's are the Pauli matrices (used also in the 
SO contributions). Note that all three components of 
the magnetic field contribute to the Zeeman term while 
only the vertical one couples with the kinetic energy 
through the vector potential. The GaAs bulk g factor 
is g* — —0.44. Finally, the Rashba and Dresselhaus SO 
Hamiltonians may be written 



hR = ( PyO x - Px&y ) , 

h D = ^L(P x a x -P y a y ) . 



(5) 
(6) 



The coupling constants \r and A^> determine the SO 
strengths and their actual values may depend on the sam- 
ple. Several experiments on quantum wells have recently 
provided valuable information about realistic ranges of 
variation for these coefficientsiii 



may use unitary transformations (as suggested in Ref. 
[T^ yielding a diagonal transformed Hamiltonian. In a 
recent worki£ we used this technique to show that the SO 
(Dresselhaus) coupling induces oscillations between up 
and down spin states when the magnetic field or the dot 
deformation are varied. Generalizing the transformations 
to consider both SO terms we define 



h = U\hUx , 
Ui = exp • 



Afi(ycr a; - xa y ) 



+ \ d (x<j x - yay) 



(7) 



Expanding in powers of the A's one finds for the trans- 
formed Hamiltonian 



2m 



+ V ext (x,y) 



N(X 2 R + \ 2 D )^ + CM\ 



g*HBB z a z 



(8) 



where we have defined the canonical angular momentum 
operator L z — xP y — yP x - Note that to 0(A 2 ), with 
A referring to both \r and A^, the Hamiltonian of Eq. 
© is diagonal in spin space. Nevertheless, the x and y 
spatial degrees of freedom are still coupled through the 
vector potential in the kinetic energy and in L z . 

A second transformation for each spin subspace of Eq. 
© may be used to obtain spatially decoupled oscillators. 
Namely, introducing a renormalized cyclotron frequency 



eB z 



+ (a!> - A 2 ) 



2m* 



(9) 



where s v — ±1 for rj =T,4, in Eqs. (5) of Ref. [13| one 
obtains the masses Mk v and frequencies Slk-q of the two 
(k = 1,2) decoupled oscillators for each spin. Anal- 
ogously, Eqs. (7) of that reference yields the eigenval- 
ues £NiN 3 7i> depending on the corresponding number of 
quanta in each oscillator (Ni, N2). For completeness of 
the presentation we repeat here the expressions for the 
latter two quantities, 



kq 



7=2 H 



CI 1 



± 



4:U! 2 0J 2 



1/2 



(10) 



where the upper (lower) sign in ± corresponds to k 
1(2), and 



(hi 



B. The analytical solution 



It is possible to obtain analytical solutions when ho ^> 



^> hz and B x 



B n 



0. In this case one 
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Sr,^g*^BB z , 



(11) 



3 





f, ^ « — 

' A A 






J 

J 

J 


i ,6=0.99 

JUl os, 




[ l A °- 95 


(N U N 2 ) = (0,0) 


i A / ? =°- 93 







1 2 

Energy £ (T) 



FIG. 1: Right panel: Evolution of the single particle ener- 
gies as a function of vertical magnetic field. Each doublet 
corresponds to different spin orientations in the transformed 
frame (Sec. II. B). The SO intensity is fixed at \ 2 D — \ 2 R — 
(1.2 x 10 _9 eVcm) 2 . The level responsible for the spin-flip 
transition when N = 7 is marked with a thick line. Left panel: 
Strength of the spin-flip excitation for different deformations 
(ui x = f3uj y ; h(ui x + toy) — 12 meV). The inset characterizes 
the transitions for N = 7 and B — of the right panel, with 
a indicating the transition between Kramers conjugates. 
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FIG. 2: Upper: spin-flip energy gap for the different levels of 
Fig. 1. The Larmor energy Tiuil is also indicated. Lower: g 
factors in absolute value inferred from the upper panel results 
as \g\ = A 3f /(jj, B B). 



As a direct application of the above results we may 
write the effective g factor for precession around a vertical 
magnetic field from the difference between the up and 
down single particle energies (A„/) with fixed oscillator 
quanta N± and N2, 



1.9 1 



fJLBB 



^bB z 



fe=l,2 



(fife-f — flkl) 



.(12) 



This equation shows that in the general case the g fac- 
tor is actually a function of the electron state (through 
the quanta), the SO coupling constants and the vertical 
magnetic field B z (through the fi's). It is also worth to 
mention that since the energy gap A s f and the modulus 
of the magnetic field (B) are positive quantities, only the 
absolute value of the g factor is determined by Eq. (|T2*|l . 



C. The transition between Kramers conjugates 

When B vanishes the full Hamiltonian fulfills time- 
reversal symmetry and, according to a well known the- 
orem of quantum mechanics, in that limit a degeneracy 
should prevail (Kramers degeneracy). As shown in Fig. 
1, the single-particle energies £niN 2 v indeed merge into 
degenerate pairs at vanishing magnetic field. These pairs 
are split by the combined action of the SO and magnetic 



field contributions and for a given (iVj^JVa) one obtains 
parallel doublets when increasing B. Depending on the 
sign of \ 2 D — \ 2 R the lower member of each doublet will 
have a given spin orientation in the transformed frame; 
namely, upwards for positive sign and downwards for neg- 
ative sign. 

If the system has good angular momentum in the in- 
trinsec reference frame (L z ), as happens in a circular con- 
finement lu x — Lu y , the Kramers conjugates at B — pos- 
sess opposite angular momenta. Therefore, the spin flip 
transition between them is forbidden since the relevant 
matrix element preserves angular momentum. On the 
contrary, when the system is deformed, for instance due 
to an anisotropic confinement oj x — (3to y , the spin flip 
transition between Kramers conjugates becomes possible 
since angular momentum is no longer a 'good' quantum 
number. This key point determines qualitatively differ- 
ent spin precessional spectra. In fact, when the transi- 
tion between conjugates is forbidden there is a gap in 
the spectrum and a non-vanishing precession frequency 
at B = (the precessional offset discussed in Ref. IT^l . 
This gap vanishes if the transition between conjugates is 
allowed due to the deformation. In the left panel of Fig. 
1 we show the evolution of the precessional peaks as the 
deformation is reduced (/3 — > 1). As discussed, the tran- 
sition between Kramers conjugates (a) switches off when 
approaching the circular case. 
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D. The g factors 

The upper panel of Fig. 2 displays the spin-flip gap for 
different levels, characterized by their oscillator quanta in 
the transformed frame. The lower panel shows the corre- 
sponding g factors obtained from A s f and the modulus of 
the magnetic field using the first equality of Eq. I|12|) . As 
in Fig. 1 a SO value of X% - X% = (1.2 x 1(T 9 eVcm) 2 as 
well as a deformation of (3 — 0.9 have been assumed. We 
note that there is a strong dependence of the precessional 
properties on the electronic state, with many cases show- 
ing a dramatic deviation from the Larmor result. When 
the number of quanta is shared asymmetrically between 
the two oscillators the g factor takes very large values at 
small magnetic fields, decreasing quite abruptly with B. 
On the contrary, when iVi = N2 there is a rather flat 
-B-dependence of the g factor and lower enhancements. 
Note also that spin-flip energies below the Larmor re- 
sult are obtained for the (0, 0) state, implying a g factor 
lower than the bare value. We have checked that other 
values of the SO couplings and dot deformations do not 
lead to qualitative variations of this behaviors although, 
obviously, the numerical values are changed. 
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FIG. 3: Numerical results for the spin-flip gap when Xr — 
and Xd = 0.5 x 10 -9 eVcm. For comparison the solid line 
display the analytical results from Eq. 1121 . 



III. CASES OF NUMERICAL TREATMENT 

When the SO coupling can not be considered weak or 
when the magnetic field points in a tilted orientation, 
with respect to the z axis, the above analytical treat- 
ment does not remain valid. One must then resort to di- 
rect numerical solution of the single-particle Schrodinger 
equation 



htfi(r,r]) = enpi{r,ri) 



(13) 



As in Ref. [13i we have proceeded by discretizing in a 
uniform grid of points, finding the orbitals and energies 
{tfi(r, rj), Ei] using matrix techniques. In terms of these 
results one can directly compute the spin-flip strength 
function, 

S prec (uj) = ^2(l~fi)fj \{<p l \(J x \ip : j)\ 2 




N= 11 


.^^(2,0) 





























1 2 
B(T) 



x S(ei 



Tiuj) 



(14) 



where i and j span the whole single particle set and the 
/i's give the orbital occupations. 



A. Vertical magnetic fields 

We have checked that the numerical solution recovers 
the previously discussed analytical limit for vertical mag- 
netic fields and weak SO couplings. For instance, Fig. 3 
compares the spin-flip gaps for cases with a weak pure 
Dresselhaus coupling having N — 7 and 11 electrons. An 
excellent agreement between the numerical data and the 



FIG. 4: Same as Fig. 3 but for Xr = 1.2 x 10" 9 eVcm and 
A D =0. 



prediction of Eq. (jl L|> is found. Note that in the numer- 
ical case discontinuous jumps in the evolution of A s f as 
a function of B are obtained whenever the ground-state 
solution implies a reordering of levels in energy. Figure 4 
displays a similar result for a pure Rashba coupling, with 
a somewhat stronger intensity. Small deviations can be 
seen with the analytical result, although the agreement 
is still quite good. Our results thus indicate that the an- 
alytic treatment works rather well for SO couplings as 
large as 1.2 x 10~ 9 eVcm, which is in the range of the 
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FIG. 5: Upper: Dependence of the spin-flip gap on the tilting 
angle of the magnetic field with respect to the vertical direc- 
tion. The thick gray line shows the Larmor energy. Lower: 
Variation of the g factors in the limit B — > as a function of 
the tilting angle 6. 



experimentally achieved values. 



B. Tilted magnetic fields 

In Fig. 5 we have analyzed the dependence of the pre- 
cessional properties on the tilting angle of the magnetic 
field with respect to the z axis, zero angle meaning per- 
pendicular magnetic held and 8 = 90° parallel B to the 
plane of motion. Note that the spin-orbit interaction is 
not invariant under rotations in the x — y plane so that 
its effects depend on the particular direction of tilting. In 
practice, however, different directions lead to only subtle 
differences, whilst the strong dependence is given by the 
angle 8. For this reason we only discuss the case of tilting 
along the x-axis. We find a rather strong dependence of 
the spin-flip gap on the tilting angle, with a maximum 
deviation from the Larmor energy for perpendicular field. 
When the tilting angle is increased a smooth energy de- 
crease in the direction of the Larmor value is seen. Actu- 
ally, for parallel orientation the results are slightly below 
the Larmor line. In the lower panel of Fig. 5 the g factors 
in the limit of vanishing magnetic field are displayed. In 
correspondence with the transition energies the largest 
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FIG. 6: Experimental spin-flip energy gaps measured in Ref. 
0for a one-electron dot. The solid line is the theoretical result 
obtained using the experimentally known u>o values while the 
dotted extension is a fit (see Sec. III.C). 



deviations from the bulk value are obtained for the per- 
pendicular direction while the parallel g factor is more 
similar to the bare factor (0.44). These results can be 
understood by noting that the SO mechanism couples 
better with the B-induced currents in the perpendicular 
geometry and, therefore, a larger influence on the spin 
precession is expected in this case. 



C. A comparison with experiment 

In a recent experiment Hanson et al& have measured 
the spin splitting in a one-electron dot by means of con- 
ductivity experiments using a parallel magnetic field. It 
is our purpose here to show that the SO-induced modi- 
fications can be the source of the observed deviation of 
the spin- flip energy with respect to the Larmor result. As 
stated in the previous section, when the magnetic field 
is alligned parallel to the plane of electronic motion the 
spin splitting recovers a Zeeman-like behavior with an ef- 
fective g-factor slightly smaller than the bulk value. This 
reduction of the spin splitting is enhanced as the spacing 
of the orbital levels is reduced, i.e., spin-orbit interaction 
induces a compression of the spin levels as w x and ui v 
become smaller. 

In Fig. 6 we display the results obtained for a cir- 
cular 1-electron dot (deformation has no significant in- 
fluence on the spin-splitting of the lowest energy state) 
with feasible values of SO coupling. Namely, we assumed 
= 0.35 x 10~ 9 eVcm, in the range of experimental 
values for GaAsp£ and X D = 0.8 x 10~ 9 eVcm. This 
latter parameter is obtained by assuming a 2DEG width 
zq ~ 60 A in the formula Ad = 7(1"/ zq) 2 , where 7 = 27.5 
eVA 3 is the GaAs specific constant w& We still need to 
input the external confinement frequency Tiwq before the 
calculation can be performed. Using for this parameter 
the measured values of the orbital level spacing, lying be- 
tween 0.96 and 1.1 meV for the range from B = to 8 
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FIG. 7: Spin precessional energies within TDLSDA. The mag- 
netic field points in the vertical (z) direction. 

Tji& one obtains the solid line of Fig. 6. For higher values 
of the magnetic field experimental values are not available 
and we have inferred the Tilhq values in order to fit the 
measured spin splittings. By assuming Hujq ps 0.5 — 0.6 
meV we obtain the dashed line in Fig. 6. Overall, the 
agreement with the measurements is rather good and, 
though this is a certainly a simplified model, we believe 
it indicates that SO coupling plays an important role in 
explaining the measured spin gaps in this system. 

D. Addition of Coulomb interactions 

The above sections have dealt with the SO-induced 
modifications of the spin precession in the absence of 
Coulomb interaction between electrons. We shall now 
estimate the role of the latter by resorting to the time- 
dependent local-spin-density approximation (TDLSDA) 
for noncollinear spins. This approach was already used 
by us in Ref.^jjfor circular dots. The reader is addressed 
to that reference for more details on this formalism. Here 
we shall only mention that the integration in time of the 
TDLSDA equations allows us to monitor the spin pre- 
cession and, in particular, to extract the precessional fre- 
quencies. Since the selfconsistent parts of the mean-field 
potential are recomputed as the system evolves in time 
one is effectively taking into account dynamical inter- 
action effects. The formalism is thus equivalent to the 
random-phase approximation (well known in many body 
theory) with an effective interaction. 

Figure 7 shows for some representative cases the pre- 
cessional frequencies in TDLSDA with SO coupling and 



deformation. A vertical magnetic field has been also in- 
cluded. We note that a qualitatively similar behaviour 
is found with respect to the preceding analytical results. 
In particular, we emphasize that at small B the preces- 
sional frequency tends to vanish and that there are dis- 
continuous jumps due to level rearrangements. It can be 
seen that, for a higher deformation (smaller /3), the in- 
dependence of the precessional frequencies is smoother, 
in agreement with the analytical model. Comparing with 
the non-interacting results there is a sizeable modifica- 
tion of the evolution in the low B range. While in the 
non-interacting scheme we obtain g factors of 20 and 6.2 
for (3 — 0.9 and 0.75, respectively, when interaction is 
included these values raise to 21.9 and 6.55. 



IV. CONCLUSIONS 

In this work we have analyzed the role of the defor- 
mation in the confinement to determine, in conjunction 
with SO coupling and magnetic field, the spin preces- 
sional properties of GaAs quantum dots. At small mag- 
netic fields the deformation closes the spin-flip energy 
gap by allowing the transition between Kramers conju- 
gate states. In practice, this implies that the precessional 
frequencies of deformed systems have no offsets at B = 0. 
The associated g factors depend strongly on the quantum 
dot electronic state and on the magnetic field direction. 
By tilting B from vertical to horizontal direction one may 
tune the g factor from large values to results close to the 
bulk one. 

When the magnetic field points in the vertical direc- 
tion and the SO coupling is weak an analytical treat- 
ment, yielding the spin-flip energies and g factors is pos- 
sible. This provides relevant insights for the analysis of 
other cases that can only be addressed with numerical 
approaches. For the case of a one-electron dot in a hor- 
izontal magnetic field we have compared the results ob- 
tained with feasible values of the SO coupling constants 
with recent experiments. We believe this comparison in- 
dicates that the SO coupling plays an important role in 
explaining the measured spin gaps in this system. For 
dots containing more electrons, the role of the Coulomb 
interactions has been estimated within TDLSDA. Size- 
able modifications of the single-particle picture have been 
obtained although, qualitatively, the main features are 
preserved. 
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